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1. Introduction 

In this work, we will study the isotropisation of a minimally coupled scalar tensor theory 
with a massive scalar field and a perfect fluid in the Bianchi type / model. 
Lets give some reasons about our choice for this geometrical framework. Although our 
present Universe seems in agreement with the cosmological principle, it could be nec- 
essary to partly release it for the early times if we want to explain the formation of 
large-scale structures. Moreover, the observations of the cosmic microwave background 
by COBE show some small inhomogeneities. These observational facts lead us to assume 
a more general geometry than this of the isotropic models i.e. the Friedman-Lemaitre- 
Robertson- Walker (FLRW) ones. The simplest generalisation is to leave the isotropy 
hypothesis and to consider the Bianchi models. About the nine Bianchi models, some 
of them accept FLRW models as exact solutions. Hence, Bianchi type I models contain 
the flat FLRW ones , the Bianchi type V ones contain the open FLRW ones, and the 
Bianchi type IX, the closed ones. At the present time, despite some more and more 
powerful observational tools, we do not know in which type of Universe we live. How- 
ever, Boomerang experiment favours a flat Universe |T] and the same conclusion could be 
drawn from the presence of inflationPHBl although this phenomenon has to be confirmed 
by deepest observations and could be compatible with other types of geometries[^. All 
these facts justify the interest of the Bianchi type / model. 

As a physical framework, we have chosen to study a scalar-tensor theory minimally 
coupled to a massive scalar field 0. The geometrical part of its Lagrangian writes as this 
of the General Relativity, which describes with a high precision the local dynamics of 
our Universe for the weak fields. Taking into accounts one or several scalar fields could 
be one of the key for a theory able to explain the physics of the early times. They are 
predicted by unification theories whose low energy limit could be a scalar-tensor theory. 
They also appear during dimensional reduction of Kaluza-Klein type theories. In addi- 
tion, we will also assume that is a massive field. The reason is that infiation could 
be the consequence of the presence of a cosmological constant whose currently observed 
value and particle physics predicted value differ from 120 orders of magnitude: this is 
the so called cosmological problem. To explain this huge discrepancy, a solution could 
be to consider that the cosmological constant is in fact a variable potential U repre- 
senting the coupling of with itself. Moreover, we will consider a Brans-Dicke coupling 
function uj between the scalar field and the metric. The theory thus described has been 
studied in jH). Here, we will generalise it by adding a perfect fluid. Associating a scalar 
field and a perfect fluid could be a way to explain the nature of dark matter, if the flrst 
one plays the same dynamical role as the second one as suggested by the quintessence 
or tracking modelslBl IZ|- In the quintessence model, the scalar field slowly rolling down 
its potential such that the ratio of its pressure and energy density, w, be a constant 
belonging to the range [—1,0]. One problem of the quintessence model is the cosmic 
coincidence problem: why the present energy density of the scalar field would be of the 
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same order as this of the matter energy density. One possible solution to this question 
could be to consider a special form of quintessence, called tracker model, for which w is 
time varying, and works like an attractor solution. Thus, late times cosmology would 
be independent of the early conditions. 

From an observational point of view, the standard model that seems to emerge 
today is the same as this described above with uj = and U = ate. However, this par- 
ticular one is far from being satisfactory (as instance, it does not solve the cosmological 
problem). Hence, more general theories leaving uj and U undetermined have to be stud- 
ied and some efforts are currently done to try and guess what could be some limitation 
on their forms and values. As instance, in [Sj, it is shown how from the observations 
one can determine the Lagrangian of a scalar tensor theory. In [0], it is demonstrated 
that scalar tensor theories can be compatible both with primordial nucleosynthesis and 
solar-system experiments with cosmological models very different from the FLRW ones. 
In [in], it is shown how new bounds on uj could be derived from future space gravita- 
tional wave interferometers, thus allowing to test scalar tensor gravity. All these works 
are related to observational cosmologies and aim to derive some satisfactory limits on 
UJ and U as we try to do it from a theoretical ground in the present paper. 

Mathematically we will study the isotropisation of the Universe in the same way 
as in [H], i.e. by associating the Hamiltonian formalism of Arnowitt, Deser and Misner 
(ADM)[TH [12] which allows getting first order dynamical system equations with the 
dynamical system methodsjin]. Our goal will be to determine the necessary conditions 
for the isotropisation of the theory described above and the asymptotical dynamical 
behaviours of the Universe at late times. The plan of this work is the following. In the 
second section, we establish the field equations of the Hamiltonian formalism. In the 
third section, we analyse their dynamics. In the fourth section, we discuss the physical 
meaning of our results. 

2. Field equations 

In this section, we will calculate the field equations of the minimally coupled scalar- 
tensor theory with a massive scalar field and a perfect fluid. The action is written: 



with (j) the scalar field, uj the coupling between the scalar field and the metric, U the 
potential and Lm the Lagrangian density of the matter. We will consider a perfect fluid 
with an equation of state p = (7 — l)p, p and p being respectively the pressure and the 
density of the fluid. For 7 = 0, 7 = 1 and 7 = 4/3 we get respectively the equation of 
state describing the vacuum energy, a dust and radiation fluid. As the first case can be 
assimilated to the presence of a cosmological constant already studied in [H], we will not 




R - (3/2 + cj(0))0'^0>2 _ ^ levrcU^ 
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consider it and will assume that 7 G [1, 2]. We will use the following form of the metric: 

ds"^ = -(N^ - NiN')dn^ + 2NidVLj + RlgijJJ (2) 

The cu* are the 1-forms defining the homogeneous space of the Bianchi type / model and 
the Qij are the metric functions. To derive the expression of the ADM Hamiltonian, we 
proceed in the same way as or ^1]. We rewrite the metric functions by using the 
Misner parameterisationfTK] : 

^33 = e-^^-^^+ 



The l3± functions describe the anisotropy of the Universe and Vt its isotropic part. We 
can derive the field and the Klein-Gordon equations of the Lagrangian formulation by 
varying the action with respect to the metric functions and the scalar field. Then using 
the Bianchi identities and the Klein-Gordon equation, we can find the conservation law 
for the energy impulsion of the perfect fiuid, T^^ = 0, and thus the well known relation 
between the energy density and the 3- volume V of the Universe: p = PqV^'^, po being 
an integration constant and V = (JliguY^'^ = e~^^. To find the Hamiltonian of the 
ADM formalism, we have to express the action with the variables /?_ and their 
conjugate momentum p+ and p_. Then, we vary it with respect to N which plays the 
role of a Lagrange multiplier. We get a constraint equation from which we can derive 
the Hamiltonian. The reader interested by full details of its derivation when no perfect 
fluid is present can find it in 0. Finally, we get: 

^2 a2 



H'=pI+ pi + 12^^^ + 247r2i?6e-6^?7 + 5e3(^-2)^ (3) 
3 -|- 2^ 

Pff) is the conjugate momentum of the scalar field and 5 a constant equal to (7 — l)po- 
It is a positive constant when 7 G [1,2] and the energy density of the perfect fiuid is 
positive. We will assume it is the case. From (jS]), we derive the Hamilton's equations: 



dH 12(f)^ 



dp^ (3 + 2uj)H 



(5) 



— f (6) 

dH dH p-^^JJ p3(7-2)f7 
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A dot means a derivative with respect to f2. If we compare them with these got when 
no perfect fluid is presentjH], we remark that only the equation ^ is modified. We will 
choose A^* = and then we calculate that A^ |II][I2] can be written^: 

N = ^^""Y (9) 
H 

The relation between the Vt time and the proper time t being dt = —Nd^, we deduce 
that t is a decreasing function of fl for all positive Hamiltonian. 

3. Isotropisation conditions and asymptotical behaviours 

In the first subsection, we rewrite the field equations with new normalised variables. In 
the second one, we study mathematically the system thus obtained. 

3.1. Rewriting of the field equations with normalised variables 



The equations (jMHI) form a first order system that we wish to rewrite with the following 
variables: 



X = H-^ (10) 
y = e-^^'VUH-^ (11) 
z = p^(P{3 + 2u)-^/^H-^ (12) 



They are independent each other's since the first one depends on if, the second one 
on H and and the third one on H, and p^. The forms of y and z show that the 
potentials U have to be positive and the Brans-Dicke coupling function u must be larger 
than —3/2 so that the variables be real. These are usual hypothesis in cosmology. The 
constraint (jH)) is then written: 

p'^x^ + R^y^ + Uz"^ + k'^ = 1 (13) 

where we have put to simplify the equations k'^ = 5e^^^~'^^^H~'^ = 5x"'y'^~^U^^'^^^ . The 
constants p and R are defined by = p+ + p"^ and R^ = 247r^i?Q. Rewriting the 
equations ([Zj) and we get: 

X = m^y'^x - 3/2(7 - 2)k'^x (14) 
y = y{Uz + 3i?\2 - 3) - 3/2(7 - 2)k'^y (15) 

z = y\3Rz - R^/2i) - 3/2(7 - 2)k^z (16) 

I In few words, this result can be recovered in the following way. We first rewrite the action under a 
Hamiltonian form, i.e. S = J {gij^-^ — NH — NiHi)d'^x, Hij being the conjugate momentum of the 
metric functions gij, N and Ni the lapse and shift functions, H and Hi the super Hamiltonian and 
super momentum. Then, by varying the expression thus obtained with respect to Hij, we derive an 
expression for dg^j/dt. Developing -l/2{-g)'^/'^Sg/Sn as a functions of the gtj and using the above 
mentioned expression, we find jHl- 
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with i = (f)U^U~^{3 + 2co')~^/^. From the first one, we deduce that x is a monotonic 
function of constant sign and then that no homoclinic orbit are allowed. 
In the [Appendix At we show that even if 5 < 0, i.e. if we consider negative energy 
density, the equilibrium is not compatible with the divergence of k. 

3.2. Mathematical study of the first order system equations 

In the first subsection, we examine the values of x and k allowing the isotropisation. In 
the second one, we look for the equilibrium points corresponding to a stable isotropic 
state for the Universe. 

3.2.1. Values of x and k compatible with the isotropisation Assuming that late times 
correspond to the divergence of the proper time t, an isotropic and stable state is such 
that tend toward some constants with dp±/dt 0. By using the expression for the 
lapse function, we calculate that it happens only when Q — oo. The Universe is then 
expanding. It corresponds to late times epoch when the Hamiltonian is positive and 
justify our assumption on the asymptotical value of t since there is no physical reason 
such that the diverging expansion of the Universe take place for a finite value of the 
proper time. Since in the same time, P± should tend toward some constants, we should 
also have P± — > 0, i.e. a; — > 0. Consequently, a stable isotropic state can be reached 
only in the plane x = of the phase space when fl diverges negatively. We will have to 
check if these two conditions do not excluded each others. 

What about the value of fc? From the constraint, we see that k^ < 1. Then, by 
considering the expression k"^ = 5H~'^e^^^~^^^ and the equation for HH issued from (jH]), 
we deduce that: 

• k"^ will tend to vanish when U > V^^ . 

• fc^ will tend toward a constant different from zero when x = if the Hamiltonian 
tends toward E = Hqc'^'^^^-'^'^^ . Then ^ 6Hq^ and U ocV-^. 

• fc^ will tend toward 1 for any potential such that U < V~''. 

The first case corresponds to an asymptotically dominated scalar field Universe. It has 
already been studied in 0. It will concern any potential tending toward a constant or 
diverging since the 3- volume will diverge when isotropisation arises and then V~'^ — > 0. 

The second case corresponds to a potential behaving asymptotically as the energy 
density of the perfect fluid, f/ oc p. It is different from trackers solutions which are such 
that asymptotically p^j, oc p, being the scalar field energy density: then, it is the scalar 
field energy density that mimics this of the perfect fluid and U < p. However, in both 
cases, the metric functions behave asymptotically as the Universe was fllled only with a 
perfect fluid whose density is increased by the presence of the scalar field. Hence, when 
U oc V~"' , dark matter and coincidence problems could be explained in the same way as 
trackers solutions try to do it. For this reason, we will name these solutions "trackers 
like" solutions. 
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The third case corresponds to a Universe dominated by the perfect fluid. Effectively, 
if we use a Lagrangian formulation to rewrite the field equations, we remark that when 
U < V~'^ , in the space field equations we can neglect the potential regarding the term of 
the perfect fiuid. Thus, the asymptotical solutions of the metric function will take the 
same form as these of a theory without scalar field. However, in the constraint equation, 
the scalar field is not negligible. We will see below that this case always implies that £ 
diverges, that is also possible in the previous case. 

In the following subsection, only the non asymptotically scalar field dominated 
cases that contain the trackers solutions, will be studied, the results of the scalar field 
dominated case being identical to these of 0. 

3.2.2. Isotropic stable equilibrium states for non asymptotically scalar field dominated 
Universe We find four equilibrium points in the plane x = 0, with k and respecting 
the constraint if k"^ = 1 — 37(2£^)~^. This last condition is not a fine-tuning. In fact, the 
expression for k contains the integration constant 6 which is hence determined by the 
constraint equation. As evoked above, it shows that when i ^ oo, k 1. Moreover, 
since /c^ and 6 are positive, we derive that £^ ^ [0,3/27] and thus eliminate two of the 
equilibrium points which are not real under this condition. The two remaining ones 
are then defined by {y,z) = {±{2RC,)''^ ^3^(2 — 7), (4£)~^7). They are real as long as 
7 < 2. The equilibrium states taking place in Q —00, they are sources for the Q 
time and sinks for the proper times t when the Hamiltonian, which can be assimilated 
to an energy, is positive. Linearising the equation (fT^ in the neighbourhood of the 
equilibrium points allows us to calculate x{fl) and, using the relation dt = —Ndfl, 
we get that asymptotically tends toward ta'^ ^ whatever i. It follows that if the 
potential behaves like the energy density of the perfect fluid, near the equilibrium, both 
tend toward zero as t~^. We have also checked that x{Q) when Q —00, thus 
showing the compatibility of these two limits necessary for isotropisation. 
In the next section, we discuss about physical meaning of these results, compare them 
with other papers and make some applications. 

4. Discussion 

In this paper, we have considered the isotropisation of a scalar-tensor theory minimally 
coupled to a massive scalar field with a perfect fiuid for the Bianchi type I model. The 
following discussion is divided in three parts. The first one contains the set of results of 
this work, the second one some comparisons with other papers and the third one, some 
applications concerning the most studied forms of potentials. 

We have seen that three cases can be distinguished depending on the fact that the 
potential is larger, smaller or behaves in the same way as the energy density of the 
perfect fiuid. For the first case, we recall the result we have got in (H], adapting its 
terms to the present paper: 
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Asymptotically dominated scalar field Universe 

When the potential is asymptotically larger than the energy density of the perfect fluid, a 
necessary condition for isotropisation of the minimally coupled scalar field with a mas- 
sive scalar field (p and a perfect fluid, whatever the Brans-Dicke coupling function u and 
the potential U considered, will he that (f)U(f,U^^{?) + 2^^)^^/^ tends toward a constant t 
with f <3. It arises at late times if the Hamiltonian is positive, at early times other- 
wise. If i the metric functions tend toward t^ ^ The Universe is expanding and 
will he inflationary if < 1. If £ = 0, it tends toward a De Sitter model. 

It includes all the diverging potentials or these tending toward a cosmological con- 
stant at late times as it could be the case for our present Universe. The first new result 
of this study concerns the potentials that mimic asymptotically the energy density p of 
the perfect fiuid. In this case, we can consider that their effect is equivalent to increase p. 

Isotropisation of trackers like theories: 

When the potential asymptotically hehaves like the positive energy density of the perfect 
fluid with an equation of state p = (7 — l)p and 7 > 1, necessary conditions for isotropi- 
sation will he that (j)'^U^U^'^{3 + 2uj)^^ tends toward a constant £^ larger that I7 and 
7 < 2. The isotropisation always arises at late(early) times when the Hamiltonian is 

2 — 1 

positive (negative). Then the metric functions tend toward the attractor t^^ , the Uni- 
verse is expanding, non-inflationary and the potential and density p tend toward zero as 
t-\ 

Finally, the last case concerns a theory asymptotically dominated by the matter. Then, 
the potential U is negligible regarding the energy density p: 

Asymptotically matter dominated Universe: 

When the potential is asymptotically smaller than the energy density of the perfect fluid, 
the necessary conditions for isotropisation are the same as for trackers like theories hut 
the quantity 0^[/|f/~^(3 + 2uj)~^ is always diverging. 

To complete these two last results, we need to determinate the asymptotical value of 
near the equilibrium. From it, we will be able to determine the asymptotical value of 
the potential relative to V~'^ and the constant £. Using the equation (0) and writing it 
near the equilibrium, we deduce that: 

Asymptotical hehaviour of the scalar field near the equilihrium: 

The asymptotical hehaviour of the scalar field near the equilihrium state when it does 
not dominate the Universe is this of the function (j) defined hy the differential equation 
(j) = S'jUUT^ when — >• —00. 
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Note that in 0, it has been shown that the corresponding equation for the scalar 
field when it dominates the Universe is = 2(I?U^{?> + 2uj)~^U~^ . All these results are 
independent from an initial state for the Universe except the sign of H which have to 
be initially positive so that the isotropisation take place at late times. Remark also that 
they can not be applied to a theory without potential since then the variable y and thus 
£ is not defined. 

In this second part, we compare our results related to the non-asymptotically dom- 
inated scalar field Universe with these of other papers. 

Hence, in [H] is studied what is called the "Quintessential adjustment of the cosmolog- 
ical constant". The quintessence phenomenon is considered for an FLRW model and 
leads naturally to a vanishing potential. It is what we observe here, the quintessential 
solutions being such that U < V~'^, since then necessary conditions for isotropisation 
imply that the potential tends toward zero at the most as It could thus solve the 
cosmological constant problem. 

In [d where the General Relativity with a perfect fiuid and a quintessential matter 
is studied for a fiat isotropic model, it has been demonstrated that the solving of the 
coincidence problem was not compatible with inflation. This result is here generalised 
to any isotropising Bianchi type / model whatever u and U since for infiation being 
present at late time we would need that 7 < 2/3, which is not the case for ordinary 
matter. 

In [iHI, it has been shown that scalar-tensor theory with a potential and a perfect fiuid 
can have as late time attractor the General Relativity. It corresponds to what we have 
found in this work since when necessary conditions for isotropy are respected the metric 

2 — 1 

functions asymptotically tend toward ts^ . Hence, when 7 = 1, we have a dust fiuid 
and the Universe tends toward an Einstein De-Sitter one with gij t^^^ as usually 
found when we consider this kind of matter. The same remark is valid when the equa- 
tion of state represents a radiative fiuid with 7 = 4/3. Then, the Universe tends toward 
a Tolman one with gij t^l'^ . 

In [in], Hyperextended scalar tensor theories (HST) with a potential for the Bianchi type 
/ model are studied. HST has the same action as (JH) but with a gravitational function 
depending on the scalar field|2n|. Some results of this last study are not changed by the 
presence of a perfect fiuid. Especially it has been shown that the Universe isotropises 
when Ge'^ tends toward a constant. Applied to the present paper, it gives that 7 have 
to be smaller than two which is the reality condition for the equilibrium points. 
If we consider General Relativity with only a perfect fiuid, it is known that the Bianchi 
type / model isotropises. In the present case, we observe that the presence of a scalar 
field add a necessary constraint, related to the range of value of the constant ^, such 
that isotropisation might arise. We also note that the asymptotical behaviour of the 
metric functions does not depends on I contrary to the case for which the scalar field 
dominates. 
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Last, from a mathematical point of view, the main difference between the case where 
the matter is absent (H] or does not dominate at late times and the case of trackers like 
theories or matter dominated theories comes from the reality condition that selects the 
equilibrium points. In any circumstances, the field equations written with new variables 
admit four equilibrium points that we will name Ei, E2, -E3 and E4. In the first case, 
reality condition selects the two first points and i belongs to a closed interval such that 
< 3. In the second case, the two last points are selected and i belongs to an open 
interval such that i"^ > I7. 

In what follows, we are going to study two well-known theories defined by 
a/3 + 2uj(l)~^ = and some exponential and power laws potential. They have mainly 
been considered for FLRW models and most of the results we will get will not be new. 
However they will permit us to test these of the present paper and to show that the 
formalism we use allow unifying them. To study each of these theories, we will proceed 
in four steps: 

(i) We calculate the asymptotical value of the scalar field if we assume that U < 
ovU> V-\ 

(ii) We respectively deduce the conditions such that U < V^^ or U > V~'^. 

(iii) We respectively deduce the conditions on i such that the Universe isotropise. 

(iv) We compare if needed, the two sets of conditions to check their compatibility. 

The first theory we want to study is defined by \/3 + 2a;0~^ = a/2 and U = e^'^. In 
[2T] . it is examined without a perfect fiuid and demonstrated that isotropisation arises 
for fc^ < 2. In [H] it is shown that it can not happen when fc^ > 6. If we suppose that 
the scalar field does not dominate at late times, we calculate that near the equilibrium 
behaves as ?)'^k~^VL. Thus whatever fc, the asymptotical behaviour of the potential 
is the same as this of the perfect fiuid energy density and hence the solution will be 
trackers like solutions. The necessary condition for isotropisation related to £ is then 
written k"^ > 37. The calculus of d(j)/dt shows that this derivative is of the same order 
or smaller than the potential and thus this solution is a true trackers one as usually 
considered in the literature. If now we assume that the scalar field asymptotically dom- 
inates the Universe and that we use the results of the asymptotical form of shows 
that our assumption is true only if fc^ < 37. The necessary condition for isotropisation 
given by the limit on £ is then A;^ < 6. Since 7 < 2, only the first inequality on k 
have to be taken into account. Hence, when the scalar field asymptotically dominates 
the Universe, the necessary condition for isotropisation is satisfied. To summarise, if 
fc^ > 37, necessary conditions for isotropisation of the Universe toward a tracker solu- 
tion such that t^^ ^ are respected whereas if k'^ < 37, necessary conditions are 
respected such that it be able to isotropise toward a dominated scalar field Universe 
with t^^ ^1211 [El- These results have been derived in for the FLRW models. 
However, in this last paper, a stable trackers solution corresponding to fc^ > 6 have 
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also been found. It does not exist in this work since then the isotropisation would be 
impossible. 

The second theory we wish to consider is defined by ^/3 + 2uj(f)~^ = ^J2 and U = (j)^. If 
we assume that the late times Universe is not asymptotically dominated by the scalar 
field, we find that (j) will tend toward (poc^'^'' (po being an integration constant§. If 
k < 0, i tends to vanish and isotropisation is not possible. If > 0, £ diverges and nec- 
essary conditions for isotropisation are respected. Thus we deduce that the perfect fluid 
will asymptotically dominate this solution, hence confirming our assumption. If now we 
suppose that U > , the calculus of the scalar field confirm it and, as shown in (H], 
the theory is able to isotropise toward a De-Sitter model asymptotically dominated by 
the scalar field when k < 0. These results are in accordance with these found in [SH] 
where it has been shown that for < 0, the solution is asymptotically dominated by 
the scalar field whereas when /c > 0, it is matter dominated. 

Some particular cases of minimally coupled scalar tensor theories with a massive 
scalar field and a perfect fiuid has already been studied in the literature. Here we 
have made an attempt to derive some necessary conditions such that an asymptotically 
isotropic stable state be reached by the Universe at late times whatever U and u and 
we have then studied its dynamical behaviour. Using these results, we have made two 
applications and checked their consistency with previous works. In a future paper, we 
hope to apply the mathematical methods of this work to the Hyperextended Scalar 
Tensor theory for which the gravitational function varies with the scalar field. 

Appendix A. Divergence of k 

Since we have chosen to consider some positive energy densities for the perfect fluid with 
moreover 7 G [1,2], then k'^ is positive and thus from the constraint, we deduce that the 
divergence of k is excluded. However, in what follows, we will consider that 6 < 0. For 
the constraint it is equivalent to write it as p'^x^ + Ry^ + Vlz^ — /c^ = 1 or to keep the 
same form as (fT3|l but with /c^ < 0. We will consider this last possibility. Then, k can 
diverge but we wish to show that it is not compatible with an equilibrium state. 
In a general manner, when x — 0, the plane where all the isotropic stable states are 
present as shown in subsection 13.2.11 we have the following relations: 

k^ Ry^ -Vlz^ (A.l) 

and then: 

y y{Uz + - 3 - 3/2(7 - 2)(1 - Ry^ - \2z^)) (A.2) 

i ^ ZRy^z - Rl2iy^ - 3/2(7 - 2)z(l - Ry^ - 12z^) (A.3) 
§ This constant does not appear in the previous appHcation since it is asymptotically neghgible. 
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The expressions (|A.2jl and (|A.3jl have to tend toward zero to reach equilibrium. For the 
first one, it will arise or 6iz + 3Ry^ -3-3/2{-f -2){1- Ry^ -12z^) 0. Let's 

study these two possibilities. 



Case 1: y ^ and (fCTl ^ 

Then, ^K.l^ implies that z diverges and ()A.2|) that yz'^ tends toward zero. Apply- 
ing these two limits to ()A.3|1 . we deduce that z ^ z"^ and thus diverges, preventing the 
equilibrium. This reasoning is also valid when (. diverges. 

Case 2: 6iz + 3Ry^ - 3 - 3/2(7 - 2)(1 - Ry^ - 12z^) and jXlI)-^ 

It means that: 



6£z + 3 + 3/2(7-2)(l - 12^^) (3/2i?7)"^ (A.4) 



k 



2 



-6iz + 3 + 3/2(7 - 2)(1 - 12z^)\ (3/27)"^ - 12^^ ^^^5^ 

By putting this last expression in (|A.3j) . we get an expression of i as a function of z. 
An equilibrium point can then be reached only for a finite value of z. But ()A.4|) shows 
that y will tend toward a constant. Thus, it will be the same for k which contradicts 
the fact that it becomes infinite. It is the same if i diverges. 

Consequently, a diverging value of k is not compatible with an isotropic state for the 
Universe at late times whatever the sign of S. 
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